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By the early 1960's advances in statistical physics had established the ex- 
istence of universality classes for systems with second-order phase transi- 
tions and characterized these by critical exponents which are different to 
the classical ones. There followed the discovery of (now famous) scaling 
relations between the power-law critical exponents describing second- 
order criticality. These scaling relations are of fundamental importance 
and now form a cornerstone of statistical mechanics. In certain circum- 
stances, such scaling behaviour is modified by multiplicative logarithmic 
corrections. These are also characterized by critical exponents, anal- 
ogous to the standard ones. Recently scaling relations between these 
logarithmic exponents have been established. Here, the theories asso- 
ciated with these advances are presented and expanded and the status 
of investigations into logarithmic corrections in a variety of models is 
reviewed. 



1.1. Introduction 

Phase transitions are abundant in nature. They are involved in the evolu- 
tion of the universe and in a multitude of phenomena in the physical, bio- 
logical and socio-economic sciences. First-order, temperature-driven, phase 
transitions, such as melting and evaporation, exhibit discontinuous changes 
in the internal energy through emission or absorption of a latent heat as the 
transition point is traversed. Higher-order transitions, in contrast, involve 
a continuous change in the internal energy. Unlike first-order transitions, 
they can involve divergences at the transition temperature T c . One of the 
major achievements of statistical physics is the fundamental explanation of 
such phenomena and, 150 years after their experimental discovery^ their 
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ubiquity ensures that their study remains one of the most exciting areas of 
modern physics. This review concerns such higher-order phase transitions. 

Second-order transitions are particularly common and examples in- 
clude fcrromagncts, superconductors and superfluids in three-dimensional 
condensed-matter physics, as well as the Higgs phenomenon in four- 
dimensional particle physics. Indeed, the edifice of lattice quantum gauge 
theory relies upon second-order phase transitions to achieve a continuum 
limit. The Kosterlitz-Thouless transitions of a Coulomb gas in two dimen- 
sions or in thin films of 4 He are examples of infinite-order phase transitions. 

A cornerstone in the study of phase transitions is the principal of uni- 
versality . This maintains that entire families of systems behave identically 
in the neighbourhood of a Curie or critical point. Nearby, thermodynamic 
observables and critical exponents, which characterise the transition, do not 
depend on the details of intermolecular interactions. Instead they depend 
only on the range of interactions, inherent symmetries of the Hamiltonian 
and dimensionality d of the system. Universality arises as the system de- 
velops thermal or quantum fluctuations of all sizes near the critical point, 
which wash out the details of interaction and render the system scale in- 
variant. This remarkable lack of dependency on the details allows us to 
understand real materials through simplified mathematical models which 
incorporate the same dimensions, symmetries and interaction ranges. Fur- 
thermore, systems with disparate physics can be categorised into common 
universality classes. Renormalization group (RG) theory provides a satisfy- 
ing, fundamental explanation of critical behaviour and universality. Indeed, 
this is one of the major achievements of statistical mechanics. 

Each higher-order phase transition is characterised by a set of univer- 
sal critical exponents. These exponents describe the strength of the phase 
transition in terms of power-laws. In the 1960's, before the discovery of 
the renormalization group, a set of scaling relations between these critical 
exponents was developed and this set is now well established and of foun- 
dational importance in the study of critical phenomena. Because of the 
importance of these power laws and the associated scaling relations, cir- 
cumstances where they are modified must be scrutinised and understood. 
In certain situations, these power laws are modified by multiplicative loga- 
rithmic factors, which themselves are raised to certain powers. Since 2006 a 
set of scaling relations between the powers of these logarithms was discov- 
ered,^ analogous to the conventional scaling relations between the leading 
critical exponents. This review focuses on these universal scaling relations 
for logarithmic-correction exponents. 
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For spin models on d-dimensional systems, logarithmic corrections of 
this type occur when the mean-field descriptions, valid in high dimensions, 
turn to non-trivial power laws for smaller dimensionalities due to the im- 
portance of thermal fluctuations there. The value of d which marks the 
onset of the importance of these fluctuations is known as the upper critical 
dimension d c . Another prominent example where multiplicative logarith- 
mic corrections are manifest is in the g-state Potts model in two dimensions. 
For q < 4 there is a second-order phase transition there, while for q > 4 
this transition is first order. Logarithmic corrections arise when q = 4. A 
more subtle example is the d = 2 Ising model with or without non-magnetic 
impurities. In each of these cases, exponents characterise the logarithmic 
corrections in a manner analagous to the way in which the standard critical 
exponents characterise the leading power-law scaling behaviour. 

The examples cited above concern well-defined Euclidean lattices. 
There, the notion of space dimensionality is crucial, the logarithmic cor- 
rections only arising at special values of d. In recent times, the study of 
critical phenomena has also focused on systems defined on networks or 
random graphs. There, instead of the system's dimensionality, a set of 
probabilistically-distributed coordination numbers characterised by a pa- 
rameter A is associated with the network. Logarithmic corrections arise at 
a critical value of A. 

Although logarithmic corrections to scaling are also encountered in the 
study of surface effects, tricritical points, the Casimir effect, and elsewhere 
in physics, the emphasis in this review is on bulk critical phenomena and 
second-order phase transitions in particular. In the following, the lead- 
ing, power-law scaling and the associated scaling relations are discussed in 
Sec ll.2| where the logarithmic-correction counterparts are also summarised. 
The standard derivation of the scaling relations in terms of homogeneous 
functions and the block-spin renormalization group is recalled in Sec ll.3l 
The logarithmic corrections are presented in Sec ll.4[ where relations be- 
tween them are derived is a self-consistent manner. Fisher renormalisation 
for logarithmic-correction exponents is discussed in Sec ll.5l In Sec ll.6l thc 
values of the correction exponents (together with the leading exponents) 
are given for various models exhibiting second-order phase transitions. 

1.2. Scaling Relations at Second-Order Phase Transitions 

Second-order phase transitions are characterised by a power-law divergence 
in the correlation length £ (the length scale which describes coherent be- 
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haviour of the system) . A consequence of this is the power-law behaviour 
of many other physical observablcs. Although second-order phase transi- 
tions are also manifest in fluids, partical physics and other arenas, in this 
exposition we adhere to the language of magnetism for definiteness, the 
translation to other fields being easily facilated. We have in mind, then, a 
system of spins s, located at the sites i of a (i-dimensional lattice, whose 
partition function is of the form 

Z = E e-? E ~ hM . (1.1) 

{*«} 

Here E and M represent the energy and magnetization, respectively, of 
a given configuration and the summation is over all such configurations 
accessible by the system. The parameters f3 = l/fc^T and h = (3H are the 
inverse temperature (divided by the Boltzmann constant) and the reduced 
external field (H is the absolute strength of an applied external field). 
The Hclmholtz free energy is usually defined as 

F L = -k B T\nZ L (T,H). (1.2) 

At h = 0, a system with N = L d sites has entropy given by 

3F 1 

Ns L = = k B In Z L + -Ne L: (1.3) 

where 

Ne L = = (E), (1.4) 

is the internal energy and (. . .) refer to expectation values. Since the pref- 
actor in Eq. (|1.2[) plays no essential role in what is to come, we drop it in 
the definition of the reduced free energy 

h = ~lnZ L (T,H). (1.5) 

With this set-up, the reduced internal energy and reduced entropy, being 
the first derivative of the reduced free energy, become essentially the same. 

According to a (modified) Ehrcnfcst classification scheme, the order 
of the transition is that of the first derivative of the free energy which 
displays a non-analycity in the form of a discontinuity or divergence. For 
magnetically-symmetric systems this occurs at ft, = and at the Curie 
critical temperature T c . We write the reduced temperature as 

t = ^j^, (1-6) 
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to express the distance away from criticality in a dimensionless way. The 
internal energy e and specific heat c arc now defined as the first and second 
derivatives of the free energy with respect to temperature, respectively, 
while the magnetization (of the entire system) m and the susceptibility x 
are defined as the first and second derivatives with respect to the external 
field. 

1.2.1. Leading Scaling Behaviour 

Since we are interested in behaviour near the critical point (t, h) = (0, 0), 
we write thermodynamic functions in terms of the reduced variables, viz. 
eL(t,h), CL(t,h), m,L{t,h) and XL(t,h), respectively. Here the subscripts 
indicate the size of the system. Then the leading, power-law, scaling be- 
haviour which describes the approach to criticality at a phase transition of 
second order is conventionally described as 

eoo (t,0) -t 1 "", eoo ((U) ~/i e , (1.7) 
CooM) ~t- a , Coo {0,h)~h- a °, (1.8) 
m oo (t,0)~t /3 for T<T C , TOoo(0, h) ~ h$ , (1.9) 

Xoo(<,0)~t-T, X oo(0,/l)~h*- 1 . (1.10) 

The numbers a, /3, 7, 5, e and a c introduced here are called critical expo- 
nents. 

The above thermodynamic functions - derivable from the partition func- 
tion - describe how the entire system responds to tuning the temperature 
and/or external field near the phase transition. To characterise local be- 
haviour within the system we need to introduce the correlation function 
and the correlation length. The correlation function is given by 

Goofat, h) ~ x-( d - 2+ rtg (xt u , ht- A ) , (1.11) 

and the correlation length in the even (temperature) and odd (magnetic 
field) sectors is usually written 

6»(t,0) U(0,h)~h~" B - (1-12) 

Again, 77, v, A and v c are critical exponents, with 77 being called the anoma- 
lous dimension. 

In addition to the above thermodynamic and correlation functions, one 
may consider the zeros of the partition function. These are the complex 
values of the temperature or magnetic-field parameters at which the par- 
tition function vanishes. For example, when h = 0, the partition function 
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in Eq. (|l.ip becomes essentially a polynomial in u = exp (—(3) when L is fi- 
nite. Moreover, this polynomial has real coefficients. As such, its zeros are 
strictly complex in the variable u and appear in complex-conjugate pairs. 
The zeros in the complex temperature plane at real h values are denoted 
by tj(L,h) and are called Fisher zerosP The set of such zeros becomes 
more dense as the lattice size increases. In the infinite- volume limit, Fisher 
zeros most commonly also lie on curves in the complex temperature planeP 
but may also be dispersed across 2-dimensional areas (typically this hap- 
pens when anisotropy is presenlP). In the thermodynamic limit they pinch 
the real temperature axis at the point where the phase transition occurs 
(namely at t = 0). 

Similarly, when /3 is real and fixed, the partition function in Eq. (jl.ip be- 
comes a polynomial in exp (—h). Zeros in the complex magnetic-field plane 
(for real values of the reduced temperature t) are denoted by hj(L,t), are 
both t- and i-dependent and called Lee-Yangzeros after their inventorsP 
In the infinite- volume limit the Lee- Yang zeros also usually form curves in 
the complex plane. In fact, in many circumstances the Lee- Yang theorerrP 
ensures that these zeros are purely imaginary. 

The impact of the Lee- Yang or Fisher zeros onto the real magnetic-field 
or temperature axis precipitates the phase transition and, in this sense, 
the zeros may be considered as "proto-critical" points^ - they have the 
potential to become critical points. Above the actual critical point T > T c , 
the linear locus of Lee- Yang zeros remains away from the real axis. In the 
thermodynamic (L — > oo) limit, its lowest point is called the Lee- Yang edge^ 
and denoted by ryi,(t). The edge approaches the real temperature axis as 
t reduces to its critical value t = and that approach is also characterised 
by a power law. Similarly, the lowest Fisher zero is denoted t\(h). For 
an infinite-sized system, this approaches the real axis as h vanishes.^ The 
leading scaling behaviour for the edge of the distribution of Lee- Yang zeros 
ar<E 

r YL (t)~t A , h(h)~h*. (1.13) 

Besides the critical indices listed above, one is often interested in the so- 
called shift exponent A s hift ■ This characterises how the pseudocritical point 
in a finite-sized system is shifted away from the critical point. The pseu- 
docritical point is the size-dependent value of the temperature t pS cudo(£) 
which marks the specific heat peak or the real part of the lowest Fisher 
zero. For a system of linear extent L it also scales as a power-law 

tpscudo(i) ~L A = hift , (1.14) 
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to leading order. We next summarise how these critical exponents are linked 
through the scaling relations. 

1.2.2. Scaling Relations for Leading Exponents 

Thus the leading power-law scaling associated with second-order phase 
transitions is fully described by 10 critical exponents a, j3, 7, 6, e, a c , 
v, v Cl A and 77 (or 11 if the exponent A s hift characterising the finite-size 
scaling (FSS) of the pseudocritical point is included). The following 8 scal- 
ing relations (9 including one for the shift exponent) are well established 
for the leading critical exponents both in the (even) thermal and (odd) 
magnetic sectors (see, e.g., Ref.[8] and references therein): 

(1.15) 
(1.16) 
(1.17) 
(1.18) 

K (1-19) 
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= 2 -a, 


2/3- 


f 7 


= 2 -a, 
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1/(2- 
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= 1, 
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= 2-<* 
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= -2 + - 






V 




v c 


= A' 
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(57 






5 - 1 



£7 

(7 + 2)(*-l) 



£7 



(1.20) 
(1.21) 



6/3 = 13 + 7. (1.22) 



The relation (|1.15j) was developed by WidorrPifll using dimensional con- 
siderations, with alternative arguments given by KadanoffP WidonP 
also showed how a logarithmic singularity can arise in the specific heat 
if a = 0, but does not have to, leaving instead a finite discontinuity (see 
also Ref. |4|10] ). Later JoscphsorP^I derived a related inequality on the basis 
of some plausible assumptions and Eq. ()1.15|) is sometimes called Joseph- 
son's relation!^ It can also be derived from the hyperscaling hypothesis, 
namely that the free energy behaves near criticality as the inverse correla- 
tion volume: /oo(i,0) ~ £,^f(t). Twice differentiating this relation recovers 
formula (| 1 . 1 5[) . For this reason, Eq. (|1.15|) is also frequently called the hy- 
perscaling relation. It is conspicuous in the set (|1.15p - (|1.22[) in that it is 
the only scaling relation involving the dimensionality d. 

The equality (j!.16[) was originally proposed by Essam and Fisher^ and 
a related inequality rigorously proved by RushbrookeP^ The relation (|1.17p 
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was advanced by Widomf^ with a related inequality being proved by Grif- 
fiths.^ Equalities (|1.16[) and (|1.17j) . sometimes called Rushbrooke's and 
Griggiths' laws, respectively^ were rederived by Abe^and SuzukP^D using 
an alternative route involving Lee- Yang zeros. Eg. 181) was derived by 
Fisher f22l with a related inequality proved in Ref.[21j. The relations (|1.19p 
and (11.201) were also derived in Refs. [18|19l22] and (|1.21[) was derived in 
Ref. [23] . The reader is also referred to Rcf.[24]. Finally the relation- 
ship between the gap exponent and the other exponents was established 
in Refs. [18119122] . 

In addition to the above scaling relations, one usually finds that 

A 8 hift = — • (1-23) 
v 

But this is not always true and in some cases, such as in the Ising model in 
two dimensions with special boundary conditions, it can deviate from this 
value A criterion for when this may happens is giverPS in Sec ll.3.4l It 
turns out that Eq. (|1.23[) may be violated when the specific-heat amplitude 
ratio is 1. 

Because of the scaling relations, only two of the exponents listed are 
actually independent. The scaling relations (|1.15[) - (|1.18[l are often listed as 
standard in textbooks, being the most frequently used. Relations (|1.19[) 
(|1.21[) involve exponents characterizing the scaling behaviour of the even 
thermodynamic functions e and c as well as of the correlation length, in 
field. Although less frequently encountered, their fundamental importance 
is similar to that of the other scaling relations. 

The formulae (|1.7[) - (|1.14[l characterise the leading behaviour of thermo- 
dynamic and correlation functions in the vacinity of a second-order phase 
transition. There are additive corrections to these scaling forms coming 
from both confluent and analytic sources. Each scaling formula is also 
associated with amplitudes, so that a more complete description of the 
susceptibility (for example) is 

X (t, 0) = r±<^ (1 + 0(t e ) + C(t)) , (1.24) 

where the amplitude r + refers to the t > symmetric phase and T_ cor- 
responds to the broken-symmetry T < T c sector. Amplitude terms such as 
these and additive-correction exponents are outside the remit of this review 
and the reader is referred to the literature! 13 ! 27 ! 
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1.2.3. Logarithmic Scaling Corrections 

Instead we focus on circumstances where the dominant corrections to the 
scaling forms (|1.7[) - (|1.14p are powers of logarithms, which couple multi- 
plicatively to the leading power laws as 



eoo(t,0) r 




(1.25) 


eoo(0, h) r 


- h e \\nhf, 


(1.26) 


Coo(t,0) r 


- t- Q |lni|", 


(1.27) 


Coo(0, h) - 


- /i- Qc |ln/i|" c , 


(1.28) 


moo(t,0) - 


-^|ln^ forT<T c , 


(1.29) 


moo(0, h) r 


- /i*|lnft| S , 


(1.30) 


Xoo(*,0) - 


- t- 7 |lntp, 


(1.31) 


Xoo(0,h) <- 


- fcT- 1 |Inft| , J 


(1.32) 


6»(t,o) - 




(1.33) 


C=o(0,/l) ^ 


- /i^ln/i^, 


(1.34) 


r YL (i) - 


-t A |lni| A fori>0. 


(1.35) 



In addition to these, the scaling of the correlation function at h = (which 
is the case most often considered) may be expressed as 

0) ~ x-( d - 2+ ^(\nxfD , (1-36) 

The above list of functions describe the salient features of a second-order 
phase transition, which is only manifest in the thermodynamic limit. The 
pscudocritical-point FSS is 

i pS cudo(i) ~ L^QnL)***, (1.37) 

It will turn out that the correlation length of the finite-sized system will 
play a crucial role and may also take logarithmic corrections, and for this 
reason we write 

^(0)~L(lnL)«. (1.38) 

Note that this allows for the correlation length of the system to exceed its 
actual length. For a long time this was thought not to be possible in finite- 
size scaling theoryPS However we shall see that this is an essential feature 
of systems at their upper critical dimensionality. Some implications of this 
phenomenon are discussed in Sec ll.7l 
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1.2.4. Scaling Relations for Logarithmic Exponents 

Over the past 5 years a set of universal scaling relations for the logarithmic- 
correction exponents has been developed j 2 l 29 l 3 ° l which connects the hatted 
exponents in a manner analogous to the way the standard relations (|1.15[) 
(|1.22j) relate the leading critical exponents. While the derivation of these 
relations is a theme of later sections, we gather them here for convenience. 
The scaling relations for logarithmic corrections are 



& 


( l + d(q-i>) if a = and cf> ^ 7r/4 
\ d(q — z>) otherwise, 


(1.39) 




= d(q - v), 




(1.40) 


Hs-i) 


= 58-i, 




(1.41) 


V 


= 7- v{2 - rf), 




(1.42) 


c 


(7 + 1)^-7) 

£ + 7 




(1.43) 


a c 


(7 + 2)(4-7) 
/3 + 7 


f 7, 


(1.44) 


8 


= dq - dv c , 




(1.45) 


A 


= P-% 




(1.46) 


Ashift 


v-q 

V 




(1.47) 



In the first of these, <fi refers to the angle at which the complex-temperature 
zeros impact onto the real axis. If a = 0, and if this impact angle is any 
value other than 7r/4, an extra logarithm arises in the specific heat. This is 
expected to happen in d = 2 dimensions, but not in d = 4, where (f> = 7r/4P 
The static scaling relations (fl~40| . (fTHjt . (ji~43j) . (fl~44| . (fi~46|) . and (fTT47|) 
can be deduced from the Widom scaling hypothesis that the Hclmholtz 
free energy is a homogeneous functiorP but the others require more careful 
deliberations, as we shall see. 

1.3. Standard Derivation of Leading Scaling Relations 

In this section we show how the leading-power-law scaling relations (| 1 . 1 5 [) 
(|1.18[) are derived using the Widom scaling hypothesis and Kadanoff's 
block-spin approach. The remaining leading-power-law relations p. 191) 
(|1.21[) are derived in a similar manner. The presentation here is necessar- 
ily elementary and the reader is refered to the standard literature (e.g., 
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Ref.[27]) for more in-depth treatments. The derivation of relation (|1.22[) 
is reserved for the subsequent section where the emphasis is on partition 
function zeros. 



1.3.1. Static Scaling 

The Widom (or static) scaling hypothesis is that the free energy / (or at 
least its singular part, which is responsible for divergences in the thermody- 
namic functions at the critical point) is a generalized homogeneous function 
(see Appendix [^) . It is convenient to express Widom homogeneity as 

f{t,h)=b- d f{V%b^h), (1.48) 

where b is a dimensionless rescaling parameter. In the renormalization 
group (RG) context, t and h are called linear scaling fields and yt and yh 
are RG eigenvalues. Differentiating Eq. (jl.48p with respect to h gives the 
magnetization and susceptibility as 

m(i, h) = b Vh - d m(b v % b Vh h) and x (t, h ) = b 2yh - d m(b v H, b Vh h), (1.49) 

while appropriate differentiation with respect to t gives the internal energy 
(entropy) and specific heat, 

e(t,h) = b yt - d f{b Vt t,b yh h) and c(t,h) = b 2vt ~ d f(b Vt t,b Vh h). (1.50) 

At h = 0, the choice b = t~ r ' Vt recovers the first expressions in each of 
Eqs. ([T7f |l - (fTTT0j) provided 

ggf - d d - y h 2y h - d 
a = , p = and 7= . (1-51) 

yt yt yt 

On the other hand, at t = 0, the choice b = h~ 1 ' Vh recovers the remaining 
parts of Eqs.([i77 | -(|LT0 l) if 

'- d ~ yt and a c = ^^. (1.52) 



d - yh yh yh 

The scaling hypothesis therefore allows one to express the six static critical 
exponents a, /?, 7, (5, e and a c in terms of yt/d and yujd. In particular, 

yt 1 Vh 



d /3(6 + 1)' d 5 + 1' 



(1.53) 



These can now be eliminated from the remaining equations (|1.51[) and (|1.52[) 
to give scaling relations (fTTT6|) . (fTTT7j) . (|TT9)) and (fOOl) . 
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The hyperscaling scaling relation (|1 . 15|) and Fisher's relation (|1.18p are 
of rather a different status than the others in that they involve the expo- 
nents v and rj, which are associated with local rather than global properties 
of the system. To derive these, we need more than the Widom homogeneous 
scaling hypothesis. 



1.3.2. Renormalization Group 

To illustrate the renormalization group, we have in mind N = L d Ising 
spins, for example, on a d-dimensional lattice with spacing a fFig |l.l[) . The 
Hamiltonian for the system is 



(1.54) 



(i,3) 



where Sj € {±1} is an Ising spin at the ith lattice site and where the inter- 
action is between nearest neighbours. In anticipation of the renormaliza- 
tion group (RG), we generalise this by introducing other locally interacting 
terms such as sums over plaquettes, next-nearest-neighbour interactions, 
etc., and we consider a reduced Hamiltonian 

H = (3H = -K t ^2 SiSj-Kh ^2 Si-K 3 ^ SiSjSkSi-K^ ^ s,s,-. . . , 

i (i,j,k,l) «j.j> 

(1.55) 

where K t = /3J, Kh = and the remaining K n are similar coupling 
strengths. Here, all of the additional terms have to be symetrical under Si — > 
Sj and <C • • • ^> indicates interactions between next-nearest neighbours. 



o o o o 



o o o o 



O O O O 



o o o o 



o o 



o o 



ba 



Fig. 1.1. In the block-spin process, the original lattice with spacing a and N spins 
becomes one of spacing ba with N/b d block spins. 



As illustrated in Fig ll.H we place the spins into blocks of length ba and 
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rescale the N' = N/b d block spins so they each have magnitude 1, 

4=.F({ Si } ie i) = ^5>, (1-56) 

iei 

where s'j labels an Ising block spin at site / of the new lattice. Here c(b) is 
a spin dilatation factor. (If all spins in a block were fully aligned then c(b) 
would be 1.) We require that the partition function of the blocked system 
is the same as the original one: 

Z N ,(W) = Z N {U). (1.57) 

This requirement induces a transformation from the point ft = {K\ , if 2 , ■ • ■ ) 
in the parameter space to another point ft' = [K'^K^, . . . ) and we write 
ft' = Rbfl. The RG approach presumes the existence of a fixed point ft* 
which is invariant under R b . If ft = ft* + 6ft and ft' = ft* + 6 ft' are nearby 
this fixed point, we linearise the RG by a Taylor expansion, 

ft* + 6ft' = Rb{ft* + 6ft) = ft* + R' b (ft*)Sft+ . . . , (1.58) 

so that 

6ft! = R! b {ft*)6ft. (1.59) 

Let the eigenvectors and eigenvalues of R' b (ft*) be given by 

Ri(ft*)Vi = \i(b)vi. (1.60) 

Demanding that two successive applications of the RG transformation us- 
ing scale factors b\ and 62 are equivalent to a single transformation with 
scale factor 6162, we have that \i{^\b-i) = \i(bi)\i(b2), so that A, is a 
homogeneous function and (see Appendix E| 

Xi(b) = b y \ (1.61) 

Vectors ft and ft* may be expanded in terms of the eigenvectors Vf. 

6ft = Uvu 5ft' = t'A- (1-62) 

i i 

The U here are called linear scaling fields. Eq. (|1.61[) then gives the linear 
version of Eq. (|1.59p to be 

i i 

and the linear scaling fields transform under the RG as 

t[ = b y %. (1.64) 
There are three cases to distinguish: 
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• If yi > 0, the associated linear scaling field £j is called relevant - it 
increases away from the fixed point under successive applications 
of the RG. 

• If yi < 0, the scaling field decreases towards the fixed point under 
the RG and is called irrelevant. 

• If yi = 0, the field ti is termed marginal. In this critical exponents 
may be continuously dependent on the parameters of the Hamilto- 
nian or logarithmic corrections may arise. 

Since the partition functions (|1.57[) of the original and renormalized 
systems are the same, the free energies near the fixed point are related by 

N'f NI {U')=Nf N {fL), (1.65) 

or, since N' = N/b d , 

f{h,t 2 M, ■■■) = ■ • ■ ) = b- d ]{W>H u V>n 2 ^H, ...). 

(1.66) 

This is precisely the Widom scaling form (|1.48j) leading to the static scaling 
relations. 

To investigate local properties, consider the correlation function of the 
block spins, 



G{tjj) = (s' iS 'j) = \ £ sWje-^, (1.67) 
{*'} 



where rjj is the distance between block / and block J in units of 
the block-spin distance. The sum is over the block-spin configurations. 
Now, for a given s'j value at a block-spin-sitc /, there is a multitude 
of sub-configurations for the set {s,} with i G I. The probability that 
the 7th blocked spin takes a particular value (say s'j = 1) is then 
£ 5(s'j — l)e~ n ^ I Z . Similarly, the probability that the blocked sys- 

tern is in any particular configuration given by Ea.(fL56l) is e ~ w ^ /Z where 



e -«'(0 = ^JJ^_jr ({a . } . 6l )) e -wW. (i. 6 8) 
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Putting this expression into Eq. (|1.67p . we obtain 

Gu = {sWj) = lE^II^-^^^e-^ (1-69) 

{s} K 



1 C 2 (b) 

Z b 2d 

{s} / \j£J 

■2 



EE< s ^> ( L71 ) 



ieijeJ 

»c 2 (b)( SlS] ) = c 2 (b)G l3 , (1.72) 

provided the two blocks / and J are sufficiently far apart that (si'Sji) sw 
(siSj) for e I and E J. If 

G. y - ~ e- r «/C, (1.73) 

then since 77 j = ry/6, we have the property that 

£' = f, (1-74) 

where £' is the block-spin correlation length, and that 

«c 2 (6)G(r,/i). (1.75) 

Demanding as usual that 0(6162) = 0(61)0(62) (homogeneity) means a power 
law, which we write as 

c(b)=b d *. (1.76) 

Here is called the anomalous dimension of the field. The choice 6 = r 
gives that G(r) ~ c _2 (r) ~ r~ 2d *. Comparing this to Eq. (|1.11[ ). we obtain 

d(j> = \{d-2 + ri). (1.77) 

Eq. (|1.64[) with ij = i and j/j = yt gives i' = b Vt t. Alternatively, ij = h 
and yi = yt gives h' = 6 yh /i. From Eq. (|1.74[) , we may write these as t£ Vt = 
t'C Vt or = h'£f Vh . Therefore we identify £ - t' 1 ^ or £ ~ /i" 1 /^, 

depending on whether h = or t = 0. From Eq. ([1.12[ ). these give 

v = — and v c = — . (1-78) 

Vt Vh 

Now, y t is related to the critical exponent a through Eq. (|1.5ip . Combining 
this with Eq. (|1.78p we obtain the hyperscaling relation (|1 . 15[) . Eq. (|1.53[) 
then leads to Eq. (|1.21[ ) for the exponent v c . 
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In summary, the scaling relations (|1.16[) . (|1.17[) t|l . 19[) and (|1.20p may 
be deduced from the Widom scaling hypothesis,^ and Eqs. (|1.15p and (|1.21j) 
can be derived from the KadanofF block-spin constructiorf^ and ultimately 
from Wilson's RGPNext we come to Fisher's scaling relation (|1.18j) . There 
are at least two ways to do this in the literature, one more common than 
the other. We describe these in turn. 

1.3.3. Fisher's Scaling Relation 

With h — 0, fixing the remaining argument of the correlation function 
p. lip , we may write 



(*) 



(1.79) 



Following the original approach used by FisherJ^ and writing the suscep- 
tibility as 

r£oo(t,o) 

Xoc{t,0)= / dtxGoofat), (1.80) 

Jo 

one obtains 

Xoo(t,0)~eoo(t,0) 2 -". (1.81) 

The leading power-laws in Eqs. (|1.10j) and (|1.12p then deliver the scaling 
relation ([TT5]) . 

An alternative argument presented in the literature^! (see also Ref. [231 
153] ) attempts to relate the correlation function to the magnetization. If the 
spins dccorrelate in the large-distance limit, one may expect that 

9oo(x,t) = (s(0)s») (s(0))(s») = n£,(t,0) (1.82) 

there. Using Eqs. ([l.llj) and (|1.12p for the left hand side and matching 
with Eq. (|1.9p on the right, one again obtains the standard scaling relation 
(|1.18p from the leading exponents. Although this approach delivers the 
correct result in this instance, we will later see that the technique delivers 
a different result at the logarithmic level in general. 

1.3.4. The Shift Exponent 

The scaling relation for the shift exponent was given in Eq. (|1.23p as the 
inverse of the correlation-length exponent v. This is immediately derived 
from Eq. (|1.48p which, with b = L and in vanishing field, implies 

c L (t,0) = L- d c{L»t,Q), (1.83) 
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having used Eq. (|1.78p . Setting the temperature derivative to zero (to 
maximise the specific heat) gives that the specific-heat peak scales as in 
Eq. (|1.14p with A s hift given by Eq. ([1.23[ ). However, it was also pointed out 
in Sec ll.2.2l that this relation does not always hold. One may ask whether 
there is a criterion for which to decide on the validity of Eq. (|1.23[) . Such a 
criterion was recently derived.^ 

Although there is no FSS theory for the impact angle of Fisher zeros, for 
sufficiently large lattice size L one may expect that <f> is approximated by 
the angle subtended by the lowest-lying zero on the real axis at the critical 
point, 

v Re[ii(L)] L-W ^ 1 ' 

Here we have used the FSS result that 7 Im[ii(L)] ~ L" 1 /". The angle is 
given by L — > oo. If, as in most cases, A s hift = then tan0 is a finite 
value. If A s hift < 1/v then tan0 = 0. But it is impossible for A s hift to 
be less than \/v - otherwise v becomes 1/A s hift- The only alternative to 
Ashift = I/* 7 , then, is A s hift > In this case tanr/> diverges as L — > oo, 
so that <p — 7r/2 and the zeros impact onto the real axis vertically. Vertical 
impact implies a symmetry between the low- and high-temperature phases, 
which in turn implies that the specific heat amplitudes on either side of the 
critical point must coincide. In other words, only in circumstances where 
the specific-heat amplitude ratio (a universal quantity) is unity is it allowed 
to violate Eq. ([L23]) . 

In the Ising model in d = 2 dimensions, vertical impact of the Fisher 
zeros and the coincidence of the specific-heat amplitudes is guaranteed by 
self-duality. In this case, Ferdinand and Fisher found A s hift = 1/v = 1 
for the square lattice Ising model with periodic boundaries,^ but other 
two-dimensional lattices with different topologies have different values of 
Ashift P 

1.4. Logarithmic Corrections 

Logarithmic corrections are characteristic of marginal scenarios (see, e.g., 
Ref.[3"5] and references therein). The hyperscaling hypothesis /co(t, 0) ~ 
£^f(t,0) fails at and above the upper critical dimension d c . While the 
leading scaling relation (|1.15p holds at d c itself, it fails for d > d c , where 
mean-field behaviour prevails. This mean-field behaviour holds indepen- 
dent of d (provided d > d c ), so the critical exponents and scaling relations 
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should also be d-independent there. At d c itself, multiplicative logarithmic 
corrections to scaling appear. 

The two-dimensional q-state Potts model has a first-order phase tran- 
sition for q > 4 and a second-order one when q < 4. The borderline q = 4 
case manifests multiplicative logarithmic corrections to scaling. The q = 2 
version of the Potts model is the Ising model. In two dimensions this has a 
logarithmic divergence in the specific heat. Although no other thermody- 
namic quantity manifests such a logarithm in this model, this scenario also 
has to be accounted for in any theory of logarithmic corrections. 

Staying in two dimensions, the Ising model with uncorrelated, quenched, 
random-site or bond disorder is another example where logarithmic correc- 
tions appear at a demarcation point. The Harris criterion^ tells us that 
when the critical exponent a of a pure system is positive, random quenched 
disorder is relevant.^ This means that critical exponents may change as 
disorder is added to the system. If a is negative in the pure system, the 
critical behaviour is not changed by such disorder. In the marginal a = 0, 
no Harris prediction is possible, and there logarithmic corrections to the 
pure model may ensue. 

Because of the ubiquity of these logarithms in critical phenomena, it 
is reasonable to seek scaling relations for their exponents in analogy to 
Eqs. (|1.15j) - (|1.22|) and (|1.23j) above. These are the scaling relations of 
Sec ll.2.3l which have only recently been developed! 2 * 29 * 30 ^ Here, these theo- 
retical developments are brought together and summarised, and their con- 
sequences are confronted with the literature. In many cases the values of 
logarithmic corrections derived in the literature using a multitude of dis- 
parate techniques are upheld. A few cases which conflict with the literature 
are highlighted as requiring further investigations. Finally, holes in the lit- 
erature are filled, pointing the way for further research endeavours into the 
future. 

The scaling theory presented in this section is entirely based on self- 
consistencies - independent of RG. The theory does not predict the ex- 
istence of logarithmic corrections; rather, when they are known by other 
methods to be present, the theory restricts their form through scaling rela- 
tions. For ab inito model-specific theories, the RG is more appropriate and 
the reader is referred to the literature! 31 * 37 tf 32E] Wegners analysis,^ in par- 
ticular, uncovered the role of marginal variables and nonlinear scaling fields. 
These were further developed by Huse and FisherP3 Excellent reviews on 
the origins of logarithmic corrections in the paradigmatic four-state Potts 
case are contained in Ref . |40) . 
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1.4.1. Static Correction Exponents 

We denote the jth Lee- Yang zero for a system of size L by 

hj(t,L) =r j {t,L)eKp(i<i>(r j (t,L))). (1.85) 

Here hj(t,L) is complex and rj(t,L) is real. The latter paramcterises the 
position of the zero along the locus of all zeros which is given by the function 
4>{r). If the Lee- Yang theorem holdsPthe angle <f>(r) is 7r/2 and the zeros 
are on the imaginary field axis. In fact the validity of the Lee- Yang theorem 
is not required or assumed in what follows (and it does not hold for the 
Potts model, for example). Expressed in terms of these zeros, the finite-size 
partition function is 

Z L (t,h) oc l[(h hj(t,L)), (1.86) 

3 

where the product is over all zeros and constant of proportionality (which 
is not displayed) contributes only to the regular part of the free energy. 
The logarithm in the free energy converts the product into a sum and 

fL(t,h) = jj52bx(h-hj{t,L)), (1.87) 

j 

up to an additive constant (not displayed). Defining 

9L(r,t) = ^^S(r- rj (t,L)), (1-88) 

3 

we may write 

f L (t,h) = [hx(h-h(r,t))g L (r,t)dr, (1.89) 



where the integral is over the locus of zeros. Since these occur in complex 
conjugate pairs, and since the density of zeros vanishes up to r = ryL,{t), 
we express the free energy in the thermodynamic limit as 

f oa (t,h)=2Re [ ln{h-h{r,t)) goo (r,t)dr. (1.90) 

Jr YL (t) 

We have assumed that critical behaviour is dominated by the Lee- Yang 
zeros closest to the critical point and that the locus of these zeros can be 
approximated by 4>{r, t) = r/>, which is a constant. We have also inserted an 
integral cutoff R. 
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The susceptibility is the second derivative of the free energy with respect 
to the external field. It is convenient to substitute r = xryh(t), so that at 
h = it is 

XcofcO) = gc,^.*) ^ (L91) 

Expanding Eq. (|1.9ip about ryi,(t)/R = 0, which is reasonable near criti- 
cality, one finds 



9oo{r,t) =Xoc(*,0)r YL (t)$ , (1.92) 

up to additive corrections in ryL (i) / R and where $ is an undetermined 
function. The ratio ryj^(t)/R is assumed small enough near criticality to 
drop additive corrections. Analogous deliberations for the magnetization 
in field give 

moo(t, h) = X oo(t, O)r YL (t)* ( — ^— ) , (1.93) 



rvL(t) 



where 



*0 ( — Kt) = 2Rc , , ^}\ ] rrdx. (1.94) 

Now, allowing /i — > in Eq. (|1.93p . and comparing to the assumed scaling 
form of Eq. (|1.29[) . one recovers the leading scaling relation (| 1 . 22[) . One also 
arrives at the first scaling relation for logarithmic corrections, namely 

A=j8-7. (1.95) 

Furthermore, fixing the argument of the function ^ in Eq. (|1.93p leads to 
t ~ h 1/A ln/i|~ A / A using Ea. ([Oo]l . so that Eq. ([Q3j) may be expressed as 

moo (t,h) ~ft 1 -*|lnft|^^f— ^-r ). (1.96) 

Next taking the limit t — > and comparing with the assumed form 
Eq. (|1.30p . recovers the leading edge behaviour (|1 .22|) . It also delivers the 
correction relation A = 8(8 — A f)/(8— 1). The former recovers the leading 
scaling relation (|1.17[) , while the latter, together with Eq. (|1.95p . gives its 
logarithmic counterpart, 

f}(5 - 1) = 88 - j. (1.97) 

Note that the assumption of logarithmic corrections to the zeros in Eq. (|1.35[) 
necessarily leads to logarithms in the other observables. If one tries to 
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omit them in other quantities (or vice versa), contradictions ensue. E.g., 
attempting to force 7 = still necessitates a logarithmic correction to the 
magnetization through a non-zero value of A in Eq. (|1.96|) . Alternatively, 
attempting to force f3 = necessitates a non- vanishing 7. By the same 
token, omission of logarithmic correction factors is the same as setting the 
relevant hatted exponents to zero and leads to contradictions. Therefore 
allowing for logarithmic corrections from the start reflects the philosophy 
of this exposition, which is based upon self-consistencies, rather than the 
Wegner approac h 1 37 * 38 ! which predicts their existence. 

We next introduce the cumulative distribution function of zeros: 

Goo (r,t) = f 9 oo(s,t)ds (1.98) 

„2 r 



= Xoo(t,0H h (t)l\—-^r\, (1.99) 

where I(y) = f, <&(z)dz. Applying integration by parts to the free energy 
(|1.90j) then gives its singular part to be 



Ut,h) = -2Re[ R G f t)ap ff . 

Jr n (t) /i-rexp(i0) 



(1.100) 



Again substituting r = xryh(t), 

/oo(t,/i)=Xoo(t,0)r YL (t)J^ » (1-101) 

where 

T,(y) = -2Re [°° ^ , (1.102) 

and we have taken the limit R/YL(t) — > 00. The internal energy and 
specific heat are given by the first and second derivatives of the free energy 
with respect to t, respectively. These are 

eoo(t,h) = XooMHlW^ 1 -^ ( — ^777 J , (1-103) 

Coo(t,h) = Xoo(t,0H L (t)tr 2 ^ C^jtA • (1-104) 
Inverting Eq. (|1.35|) . we may write 

t~4 L \\nr YL \-%, (1.105) 
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which with Eq. (|1.31[) gives 



Xoo(*,0) ~r Y *|lnr YL |^+r (1.106) 
Together, these give 

eoo (t,h) = 4-*"*|lnr YL | i2 ^+^ (—) , (1.107) 



Coo(^^) = 4l A A |lnr YL |^^+^ ( — j . (1.108) 

Comparing Eq. (|1.108|) with the form (fL~2"7]l , one obtains a = 2 + 7 - 2 A 
and a = 7 + 2 A. From Eqs. (|1.17|) and ()1.22|) . the former is the standard 
scaling law (|1.16[) . From Eq. (|1.95[) . the latter can be expressed as a third 
relation between the logarithmic-correction exponents, namely 

a = 2/3 -7. (1.109) 

Eqs. (|1.107|) and (|1.108j) may be rewritten as 

eoo^fc) = h'-i-^lnhl^^ (^) , (1.110) 

Coo (t,h) = / l 2 -*-i|ln/ i | i2 ^+^ (—) . (1.111) 

Letting t — > so that tyl(0 — > 0, one finds 

eoo(h) = h 2 -Z-i\lnh\ td £ A+ * / , (1.112) 

Coo W = /i^i-illn/il 12 ^^. (1.113) 

From the leading behaviour one recovers Eqs. (|1.19[) and (|1.20[) . From the 
logarithmic corrections, one obtains the counterpart scaling relations ()1.43j) 
and ([PI)l . 

1.4.2. Hyperscaling for Logarithms 

The logarithmic analogue of the hyperscaling relation (|1.15|) has a rather 
different status than the other critical exponents, in that it necessitates 
consideration of finite-size effects (this aspect will be further discussed be- 
low). Consider, therefore, a system of finite volume N = L d . The finite-size 
scaling (FSS) of Lee- Yang edge is given by 

ri (L) (1 . 114) 



r YL (i) VCoo(t,0) 
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where £l(0) is the correlation length of the finite-size system at t = 0. To 
allow for logarithmic corrections to this quantity too, we assume the form 
(|1.38[) . For finite systems, the cumulative density of zeros is the fractional 
total of zeros up to a given point 

G L ( rj (L)) = ^^. (1.115) 

Fixing the ratio r/ryi,(t) in (|1.99|) . and using the scaling relations previ- 
ously derived, and then taking the limit t — > 0, one arrives at an expression 
for the critical cumulative distribution function, 

GooM) ~ r ^r\\nr\ & ~^ A . (1.116) 

At t = and for sufficiently large L, the FSS expression f| 1 . 1 1 5[) must 
converge to the infinite- volume expression (jl.1161) . Equating them gives 
the FSS of the first Lee- Yang zero at criticality to be 

n(L) ~L-^ (lnL) A ~^ . (1.117) 

Inserting (fL~55]l , (TOS]) . (fT35|) and (|1.117p into (|1.114[) recovers ([TJ5j) and 
yields the logarithmic equivalent to the hyperscaling relation, namely 

a = dq-di>, (1.118) 

The logarithmic scaling relations (|1.97[) and (|1.109l) but not f| 1 . 1 1 8[) can 
be derived starting with a suitable modification to the phenomonological 
Widom ansatz,E2E31 namely 

f(t,h) = b- d f(b v < ' \\nbfn,b yh (lnb) ih h). (1.119) 

Following the approach of Sec ll.3.11 this recovers the static scaling equa- 
tions (|1.16l) . (|1.17|) . (|1.19[) and (|1.20l) as well as their logarithmic counter- 
parts (fTlOj) . (|L4Tj) . (fT43|) and ([L44jl with 

/3 <5<5 2/3- 7 , 

* = f " flTT*) = (L120) 

and 

w = TTi = W (L121) 

The logarithmic counterparts of Eqs. (|1.51[) and (|1.52[) are then 



" = A P = Pyt + Vh, l=-iyt + 2y h (1.122) 

yt 
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and 

5 = d—, e = y t + ey h , a c = 2y t - a c y h , (1.123) 
Vh 

while 

A = Ay t - y h . (1.124) 

1.4.3. Logarithmic Counterpart to Fisher's Relation 

With multiplicative logarithmic corrections, the correlation function is 
given by Eq. (|1.36[) . Fixing the argument of the function D, we may rewrite 
this as 

Qoo{x,t,o) ^ u(t,or id - 2+ri) (iHo.it, o)rD {j^roj) ■ (L125) 

Following Fishei^ and writing the susceptibility in terms of the correlation 
function (Eq. (|1.80l) '). one obtains 

Xoo(t,0) ~ £oo(t,0) 2 - ,, (ln{; oo (t,0)) f >. (1.126) 

As in Sec ll.3.3l the leading power- laws in Eqs. (|1.31j> and (|1.33j) then deliver 
the scaling relation (|1.18|) . Matching the logarithmic corrections exponents 
too yields 

f) =-7 — 1/(2 — 7?). (1.127) 

This approach was used in Ref.[32] for the case of the d — 2 four- 
state Potts model. An alternative argument presented there (see Sec ll.3.3]) 
assumes that the spins decorrelate in the large-distance limit, 

Goo(x,t) = (,s(0) S ») ^ (s(0))( S ») = m^(t,0). (1.128) 

Using Eqs. (jl.33p and f| 1 . 1 25[) for the left hand side and matching with 
Eq. (|1.29p on the right, one again obtains the standard scaling relation (|1.18p 
from the leading exponents (see also Refs l23l33)) . However, it is interesting 
to note that matching the logarithmic exponents gives, instead of (| 1 . 1 271) . 
the relation f) = 2/3 + i>(d — 2 + rj). From the logarithmic scaling relation 
(fl~40|) this would mean 

f) = dq + *f-v(2-ri), (1.129) 

which is, in general, different to Eq. (|1.127[) . 

When q vanishes, as is the case in the d = 2, four-state Potts model P^l 
for example, this is actually identical to Eq. (|1.127[) . To decide between 
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Eqs. (|1.127[) and (|1.129[) . we need a model with non-zero q value. Models at 
the upper critical dimension provide suitable dctcrminators and one finds 
that the Eq. (|1.129[) fails in these cases. Eq. ([1.127p holds in each case we 
have examined. 

1.4.4. Corrections to the Logarithmic Scaling Relations 

Although it holds in most models which manifest multiplicative logarithmic 
corrections, there is an immediate and obvious problem with scaling rela- 
tion (|1.118j) when it is confronted with the Ising model in two dimensions, 
which has been solved exactly. There, the specific heat has a logarithmic 
divergence so that a — and a = 1. There are, however, no logarithmic 
corrections in any of the other thermodynamic or correlation functions in 
this model, and, since v = q = 0, Eq. (|1.118[) fails immediately. The relation 
(II. 118P also fails in the uncorrelated, quenched, random disordered version 
of the Ising model in two dimensions, where q = 0j 45 l 46 l a = 0, v = 1/2 
ancPSMl 

Coo0,0) ~ln|lni|. (1.130) 

This famous double logarithm in the specific heat of the diluted Ising model 
in two dimensions has been the source of great controversy throughout the 
years. This controversy has only recently been convincingly resolved, partly 
with the aid of logarithmic-correction theories} 52 * 53 ^ which should be able 
to account for them. In this section, besides resolving the puzzle as to the 
value of a in the pure and random Ising models in d = 2 we will see how 
this double logarithm emerges quite naturally from the general scheme. 

Since the problem is associated with the even sector of the model 
(namely with the i-dependency of the specific heat), one may argue that 
a Lee- Yang analysis, which focuses on an odd (magnetic) scaling field, is 
not the best approach to fully extract the general relationship between the 
correction exponents appearing in Eq. (|1.118[) . One may appeal to complex- 
temperature (Fisher) zeros for further insight, as they are appropriate to 
the even sector. We will see that the puzzle is resolved after consideration of 
two special properties of the pure and random two-dimensional Ising mod- 
els, namely the vanishing of a and the angle at which their Fisher zeros 
impact onto the real axis. 

An FSS theory for Fisher zeros is obtainecP by writing the finite-size 
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partition function in terms of the scaling ratio £t / £ 



^M) = Q(-%^). (1-131) 



,foo(*,0). 

This vanishes at a Fisher zero. Labeling the j th finite-size Fisher zero as 
tj(L), one therefore has 

CTT°' 1(0) ' (L132) 

where QJ 1 ^) is the j th root of the function Q. Using EdS- lfOS) and (|T3"g]) , 
this equation can be solved to give the FSS form of the jth zero, 

tj (L) ~ L~» (In L) ^ . (1.133) 

So far, no assumptions other than the validity of FSS have been involved. 
Note that this is the same scaling form as Eq. (|1.47[) . which is sensible, since 
the real part of the lowest zero is also a pseudocritical point. 

The full expression for the scaling of the j th Fisher zero i j 7 l 42 l 54 l a func- 
tion of a fraction of the total number of zeros Af. Since this is proportional 
to the lattice volume N = L d , Eq. (|1.133|) is more appropriately written as 

1 O — Q 

tj (L) ~ fc^) ^ (in (^r)) " exp^r^Z))), (1.134) 

where (f>j(rj(L)) is the argument of the jth Fisher zero. 

In all cases known from the literature, the Fisher zeros for isotropic mod- 
els on homopolygonal lattices fall on curves in the complex plane. They 
impact onto the real axis along so-called singular lines lilSSEH We a i so as _ 
sume this scenario in the present case, and we denote the impact angle onto 
the real axis in the thermodynamic limit by r/>. 

Similar to the Lee- Yang case, we may write the finite-size partition 
function in terms of Fisher zeros, 

Z L (t, 0) oc H (t - tj{L)) (t - t*(L)) , (1.135) 
i=i 

where we have been careful to identify tj(L) and t*{L) as complex conjugate 
pairs. Provided that the M. ex M zeros which dominate scaling behavior in 
the vicinity of the critical point are described by the scaling form (|1.134|) . 
and differentiating appropriately, one finds the FSS for the specific heat at 
t = to be 

M 

C L (0) - -L- d Rc^i7 2 (L). (1.136) 
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We wish to compare this expression with that resulting from a direct 
application of the modified FSS hypothesis (|1.114j) to the specific heat, 
which yields 

c L (o) WMW). (L137) 



Coo(t,0) \Zoc(t,0) 

Fixing the scaling ratio on the right hand side so that t ~ 
L- 1 / l/ (ln£)( p - < ?)/ l/ , one obtains from (jl~27)) the FSS behaviour 

C* L (0) ~ L»(kiL) & - a ^. (1.138) 

We now match Eq. p.l36j) to Eq. ([1.138|) . 

In the case where a ^ 0, Eq. (|1.136p gives 

C L (0) - L» {InLy 2 ^ , (1.139) 

and comparison with Eq. (|1.138|) leads to the recovery of the previously 
derived logarithmic scaling relation (|1,118[) , 

If, however, a = 0, the FSS expression (|1.136[) for the specific heat 
becomes 

^cos(2^(L)) ( j-1/2 



^-^-T^rV-Lt-)) (L140) 

For large L and close enough to the transition point that (f>j(rj(L)) ~ <f>, 
the cosine term becomes a non- vanishing constant if ^ 7r/4. This is 
what happens in the pure Ising model in d = 2 dimensions, where (f> = ir/2 
is assured by duality on the square lattice^ as well as for its symmetric 
random-bond counterpart P Universality of <f> ensures this happens for other 
lattice configurations too and continuity arguments lead one to also expect 
^ 7r/4 in the general random-bond and random-site two-dimensional Ising 
case. 

In these cases, application of the Euler-Maclaurin formula gives that 
the leading FSS behavior when a = is 

In the thermodyamic limit, one may simply replace L and Cl(0) by t 
and Coo(t, 0) in Eq. (|1.14ip . respectively. Then, comparing Eq. (|1.138[) with 
Eq. (|1.14ip and using standard hyperscaling (|1.15[) , one finds 



a = 1 + dq — dv. 



(1.142) 
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This expression replaces Eq. (|1.118p when the model has both a = and 
4> / 7r/4. Thus the logarithmic scaling relation (jl.39j) is established. Fur- 
thermore, we prefer to write Eg. ([1. 401) in terms of q and v rather than in 
terms of a as in Eq. (|1.109|) because, unlike a, q and v the q~v combination 
does not exhibit the subtleties discussed here. 

We can perform similar considerations for the Lee- Yang zeros, replacing 
Ea. (li~T3il by 

Z L (h) = Q( } L ^\, \ . (1.143) 



We find 

5 = d(q - u c ), (1.144) 

except when 7 = and the impact angle is not 7r/4, in which case an extra 
logarithm appears and Eq. (|1.144p becomes 

5 = 1 + d(q - v c ). (1.145) 

That the impact angle is 7r/2 is guaranteed by the Lee- Yang theorem, but 
we know of no cases with 7 = 0, so we omit Eq. (|1.144[) from Eq. (|1.45p . 

To summarize so far, we have explained where the standard scaling re- 
lations come from using the standard approaches of Widom scaling and 
the Kadanoff block-spin RG. These standard scaling laws for the leading 
critical exponents are well established in the literature. Analogous rela- 
tions for the logarithmic corrections (|1.39p - (|1.46p have also been derived 
using a self-consistency approach. Next we confront the logarithmic scaling 
relations with results from the literature in a variety of models. Where 
predictions for, or measurements of, these logarithmic corrections already 
exist in the literature, we can test the new scaling relations. Indeed, we 
will show that they are upheld (barring at least 3 exceptional circumstances 
which require further research). Where there are gaps in the literature re- 
garding the values of the logarithmic corrections, we can use the scaling 
relations to make predictions for them. These new predictions need to be 
independently tested in future research. 

1.4.5. The Logarithmic Shift Exponent 



Finally we mention that the finite-size scaling of the pseudocritical point 
may be determined using the Widom static hypothesis (|1.119p modified to 
include logarithms. In a similar manner to Sec ll.3.4l one finds that the 
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specific heat peaks scales as Eg. (|1.37p with 

Ashift = -Vt = V —- (1-146) 
v 

This is the final scaling relation for logarithmic corrections Eg. (II. 471) . 

1.5. Fisher Renormalization for Logarithmic Corrections 

Fisher renormalization concerns systems under some form of constraint. 
For such systems, the critical exponents take values which differ from their 
"ideal" counterparts. The systems we have been dealing with so far are 
"ideal" , in the sense that they are not subject to constraints of this type. 
Typically, the theoretical power-law divergence of the specific heat in an 
ideal system, for example, is replaced by a finite cusp in a "real" experimen- 
tal realization. Fisher^l explained this as being due to the effect of hidden 
variables and established elegant relations between the ideal exponents and 
their constrained counterparts. 

Phase transitions which exhibit Fisher renormalization include those in 
constrained magnetic and fluid systems (e.g., with fixed levels of impuri- 
ties), the superfluid A transition in 3 He- 4 He mixtures in confined films 
the order-disorder transition in compressible ammonium chloride, 1 ^^ the 
critical behaviour at nematic-smectic-A transitions in liquid-crystal mix- 
tures^ and in emulsions P^l 

When multiplicative logarithmic corrections are present in the ideal sys- 
tem, these migrate to the Fisher-renormalized system in a non-trivial man- 
ner which has recently been determined.^ Here, we summarize the Fisher 
renormalization for the logarithmic exponents. We later use this scheme 
to deduce the scaling behaviour of lattice animals at their upper critical 
dimension. 

For a system under constraint, Fisher established that if the specific-heat 
exponent for the ideal system a is positive, it is altered in the constrained 
system. The magnetization, susceptibility and correlation-length critical 
exponents are also changed. If the subscript "X" represents the real or 
constrained system, its critical exponents are related to the ideal ones by 
the transformations 

ax = z , Px = , lx = z , v x = z ■ (1.147) 

1 — a 1 — a I — a 1 — a 

The exponent 5 and the anomalous dimension r\ are not renormalized: 

S x = S, and i] X = V- (1.148) 
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Expressing Eqs. (|1.53[) in terms of the above exponents, and Fisher renoma- 
lizing, gives the RG eigenvalues for the constrained system as 

Vt x = (! - a )Vt, and y hx = y h . (1.149) 

One may then use the scaling relations (|1.2ip . (|1.22|) . (|1.23|) and (|1.52l) to 

determine the Fishcr-rcnormalization formulae 

A a a . 

Ax = - , v cX =v c , e X = f- + -r, a cX = a c -2—. (1.150) 

1 — a A A 

as well as 

A s hiftx = d — A s hift. (1.151) 
These formulae have two attractive properties: 

• If the ideal exponents obey the standard scaling relations then the 
Fisher renormalized exponents do likewise. 

• Fisher renormalization is involutory in the sense that Fisher renor- 
malization of the constrained exponents returns the ideal ones. 
This means that two successive applications of the transformation 
gives the identity. 

Fisher renormalization applied to the most commonly encountered log- 
arithmic correction exponents results irP^l 



ax 


1 


a 

3 

— a 


(1.152) 


hx 


= h- 


Pa 
1-a 


(1.153) 




= 7 + 


7a 
1-a' 


(1.154) 


t>x 


= z> + 


va 
1 — a 


(1.155) 



As for the leading indices, no Fisher renormalization occurs for the 
logarithmic-correction exponents for the in-field magnetization or for the 
correlation function, since it is defined exactly at the critical point. Simi- 
larly q is unchanged. 

5 X = S, fix =fj, q x = q- (1.156) 

Each of these obey the scaling relations for logarithmic corrections, and 
Eqs.(EIU|) and (fTT2T|) give 

Vtx = (1 - oe)y t — a, (1.157) 
Vhx = yh, (1.158) 
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while Eas. (fn5]) - (fT37p also lead to 

aA 

e x = e-a-— , (1.159) 

(1.160) 

(1.161) 
(1.162) 
(1.163) 



OtcX 


= a c 


VcX 


= Vc, 


A x 


= A 


shiftX 


= (1 




Ac 



1-a 



Fisher renormalization is also involutory at the logarithmic level: re- 
rcnormalizing the Fisher-renormalized logarithmic exponents returns their 
ideal counterpart. In other words, Fisher renormalization is its own inverse 
at the logarithmic as well as leading level. 

Fisher renormalization at the logarithmic level has been tested in lattice 
animals and the Lee- Yang problem at their upper critical dimensions, which 
also led to new predictions for logarithmic corrections are made. These are 
discussed in Sec ll.6.6l below. 

1.6. Logarithmic Correction Exponents for Various Models 

In this section we confront the logarithmic scaling relations (|1.39j) - (|1.46j) 
with the literature for various models on a case-by-case basis. 

1.6.1. q = 4 d = 2 Potts Model 

The leading critical exponents for the 4-state Potts model in d = 2 dimen- 
sions arc™ 

2 1 7 2 

a=-, P = -, 7 =g, 6 = 15, „=-, 

1 4 8 8 A 5 

V=~ A , e=-, a c =-, u e =-, A = I . (1.164) 

The hitherto-known logarithmic-correction exponents ar d 32 * 40 * 6 -^ 

a = —1, B = — , 7 = — , 5 = . z> = — , 77 = — . (1.165) 

8 4 15 2' ' 8 v ; 

FSS of the thermodynamic function d 32 * 66 ! means that 



4 = o. 



(1.166) 
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The correction relations (fT3"gj) - ([T12]) therefore hold, while (fTT^j) gives A = 
5/4 for the leading scaling of the Lee- Yang edge. We can also use the 
remaining static scaling laws (|1.43p f|1.45[) to predict 

6 = "30' ac = _ 15' " c= 30' A = ~8- (L167) 
1.6.2. O(N) c(>* Theory 

The upper critical dimension for 0(iV)-symmetric r/) 4 . theories - of crucial 
importance to the Higgs sector of the standard model - is d c = 4. Here 
hyperscaling fails and the leading critical exponents take on their mean-field 
values which are 

a = 0, (i = 7 = 1, (5 = 3, v ■ = -, 

»7 = 0, e= |, a c = 0, i/ c = i A = |. (1.168) 

The RG predictions for the corrections are already known to b d 44 * 67 * ^ 
4 - JV - 3 N + 2 - I 

s= ivT8' p = W+s' ^ = ivT8' * = s' (L169) 

iV + 2 t 1-N 1 

r/ = ' A = ^Ts' q = l- (L170) 

The correction relations (|1.39p - (|1.42p and (|1.46p therefore hold and from 
formulae (|1.44[) - (|1.45|) we predict 

10 -N . 4- AT „ 1 
£= 3(ivT8)' ac = ^T8' " c= 6- (L171) 

In the four-dimensional Ising case it is knowrJ^iMI th a t the impact angle 
for Fisher zeros is <f> = 7r/4. Assuming the same for its O(N) generaliza- 
tion, means that Eq. (|1.142p does not follow from Eq. (|1.140p and Eq. (|1.118|) 
remains valid there instead. The same is expected to be true for O(N) the- 
ories with long-range interactions, which we next address. 



1.6.3. Long-Range O(N) <t>% Theory 

The introduction of long-range interactions alters the universality class of 
O(N) spin models. If the interactions decay as x~^ d+(T \ the upper critical 
dimension becomes d c = 2a. The critical exponents for the A^-component 
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long-range system are 1 70 ! 71 ! 

a = 0, p = \, 7 =1, S = 3, v=-, (1.172) 
Z a 

V = 2-a, e=|, a c = 0, v c = ±, A = |, (1.173) 

which obey the leading scaling relations. The Privman-Fisher form for the 
free energy^ gives the RG values for the logarithmic corrections to be 

„ _ 4- N 3 „ _ iV + 2 

a ~ iV + 8' ^ ~~ + 7 ~ iV + 8' 

5=1 t)= Ar + 2 N , n = 0. (1.174) 
3' a(N + 8)' ' K ' 

The correction relations (|1.40p - (|1.42[) are obeyed, and the remaining re- 
lations ([1755)1 and ([T733")) - (fT74"g)l predict 

10-7V „ _4-iV „ _6-rr i_l-^V 
e ~ 3(jv + 8)' ac_ ivT8' ^ C ~T2V' 77T8' 9 ~27' 

(1.175) 

The prediction q = l/2er recovers the known valueP q = 1/4 for O(iV) 
r/>| theory in the a — 2 case. It also leads to agreement with long-range 
Ising FSS in two dimensions when a = lP^ The remaining predictions for 
long-range systems have yet to be verified. 

1.6.4. Spin Glasses in 6 Dimensions 

Spin glasses, percolation, the Lee- Yang edge, and lattice animals problems 
are each related to <fi 3 field theory. The leading exponents are^HZH 

o=-l, 13=1, 7 = 1, 8 = 2, i/=i, r? = 0, (1.176) 

and obey the standard scaling relations provided 

A = 2, e = l, Oo = ~, v c = - A , (1.177) 

predictions which still need numerical verification. Ruiz-Lorenzo used RG 
methods to derive the critical scaling exponents of the correlation length, 
susceptibility and specific heat for these models at their upper critical di- 
mensions as^ 

„ 2(2o-3a) 2a „ 5a . „ oN 

a = Ab-a ' 7 = 4b— a U = 6(46^)' (1 ' 1?8) 

where the values of (a, b) depend upon which problem one is considering. 
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In the m-component spin-glass case, the upper critical dimension is 
d c = 6 and values of (a, b) are (—4m, 1 — 3m). There, the leading exponents 
are given by Eq. ([1.176p while Eq. (|1.178p gives the logarithmic correction 
exponents as 

3m + 1 2m 5m 

a = ~o T' 7= 9 v v = TTn TV 1179 

2m — 1 2m — 1 6(2m — 1) 

These satisfy the scaling relations (|1.39j) - (|1.46|) provided that 

- 1 + m - 1 — 3m „ m ,11 + 2m 

' ~ 2(1 - 2m)' = 4(1 -2m)' V = 3(2m - 1) ' € = 2 1 -2m' 

1 3 + 7m „ 1 5771 — 3 - 1 + 5m 1 , . 

ac = IT^' ^= 242^1' A= 2(l-2m)' 9= 6' (L180) 

Independent numerical investigations of these logarithmic corrections are 
required. In particular, Ruiz-Lorenzo's prediction for the finite-size 
correlation- length correction exponent is q = 1/3. 75 This is the first in- 
stance where we encounter a clash with results in the literature regard- 
ing the exponent q, the other two being in the percolation and the Lee- 
Yang/lattice-animal problems. These cases are discussed next. 

1.6.5. Percolation in 6 Dimensions 

The percolation problem at its upper critical dimension of d c = 6 has^EU 
the same leading, mean- field critical exponents as in Eqs. (|1.176j) and 
(|1.177|) . With (a, b) = (—1,-2), in the percolation case the following cor- 
rection exponents are knowrP^S] 

2 ~ 2 2-2 5 1 

&=-, B=-, 7=-, 5=-, 7>=— , 77=—. (1.181) 

7 7 7 7 42 21 y 

The scaling relations for logarithmic corrections now allow the prediction 

<=\ &c=^, "c = ^, A = 0, g=i. (1.182) 

While other works^ contain an implicit assumption that q = 0, Ruiz- 
Lorenzo's prediction^ for this quantity is q = 1/3. Again, further investi- 
gations are needed. 

1.6.6. Lee- Yang Problem in 6 Dimensions and Lattice Ani- 
mals in 8 Dimensions 



The lattice- animal problenPHM] [ s closely linked to the Lee- Yang edge prob- 
lemP The former has upper critical dimensionality 8, while for the latter 
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it is 6. The mean- field values of the critical exponents for the Lee- Yang 
edge problem are again given by Eqs. (|1.176|) and (|1.177j) . while the values 
(o, b) = (—1,-1) for Lee- Yang singularities^ lead to 

a = -\, T=| *=1. (1-183) 
The scaling relations for logarithmic corrections then predict 

/3 = 0, <5=i r) = i e = 0, 

1 1-2 1 

a c = --. v, = -, A = --. <?=-. (1.184) 
3 9' 3' q 6 v ' 

Fisher-renormalizing both the mean-field leading critical exponents and the 
logarithmic corrections, one obtains (omiting the subscript X) 

V = 0, A = l, e =i ac=~ i/ c = i (1.185) 



and 



- 1 fl 1 - 1 x 1 - 1 

3' M 3' r 3' 3' 9' 

A = 0, c=i a c = i ;> c = i (1.186) 

together with 

q=\- (1-187) 

These deliver our theoretical predictions for the lattice-animal problem at 
its upper critical dimensionality d = 8. 

A lattice animal is a cluster of connected sites on a regular lattice. Vari- 
ants include clusters of connected bonds as well as weakly embedded and 
strongly embedded trees. It is believed that these models belong to the 
same universality class. The enumeration of lattice animals is a combina- 
torial problem of interest to mathematicians. In physics, they are closely 
linked to percolation and clustering in spin models. In chemistry they form 
a basis for models of branched polymers in good solvents. Lattice animals 
linked by translations are considered as belonging to the same equivalence 
class, and as such are considered to be essentially the same. One is inter- 
ested in Zn, the number of animals containing N sites, and the radius of 
gyration Rn, which is the average distance of occupied sites to the centre 
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of mass of the animal. Allowing for logarithmic corrections, these behave 
asP 

Z N ~n N N-"{lnN) d , (1.188) 
R N ~ N v QnNf. (1.189) 

The exponent 9 may be identified with 3 — a, which, from the above Fisher- 
rcnormalized values is 5/2. There is a famous scaling relation due to Parisi 
and Sourlas which predicts that 9 and v are related bjE^ 

6 = (d-2)v + l. (1.190) 

This is essentially hyperscaling (|1.15[) with the d dimensionally reduced 
by 2 and has been numerically verified in dimensions d = 2 to d = 9.^1 
Identifying 9 = a The scaling relation (| 1 . 39[) may be written 

9 = (d-2)(q-i>) = 6(q-P), (1.191) 

having reduced the dimensionality term from d = d c = 8 to 6, appropri- 
ately. This is the logarithmic counterpart to the Parisi-Sourlas relation 
([LITO]) . The above values a = 1/3, q = 1/6, z> = 1/9 satisfy Eq. (frTOT|) 
by construction. An alternative set of values in the literature^ is a = 1/3, 
q = 1/3, v = —1/72. This set comes directly from an RG-based calculation 
and does not satisfy the new scaling relation. Moreover, a recent numerical 
study though not conclusive, indicates that the set of exponents devel- 
oped here is more likely to be the correct one. Clearly more research is 
needed to explain the disparity of this set with the RG. 

1.6.7. Ising Model in 2 Dimensions 

The Ising model two dimensions has critical exponents 

1 7 

a = 0, r3=g, 7=4> <*=15, f = l, 

»7 = T. A = 0, e = 0, a c = 0, u c = 0. (1.192) 

The specific heat has a logarithmic divergence with temperature in this 
model so that a = 1. However, this well-studied model is unusual in that 
there are no other logarithmic divergences and all the remaining hatted 
exponents vanish. The subtle nature of the emergence of this unusual be- 
haviour was explained in Sec ll.4.41 a logarithm arises from summing over 
the Fisher-zero indices, a summation which is necessary because these zeros 
approach the real-temperature axis at an angle other than 7r/4. In fact the 
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angle of approach is ir/4. This is because the self-dual nature of the two- 
dimensional Ising model assures a symmetry between the high- and low- 
temperature sectors, which demands that the impact angle </> = ir/2. 

1.6.8. Quenched- Disordered Ising Model in 2 Dimensions 

The uncorrelated, quenched, random removal of sites or the randomisa- 
tion of bond strengths on a lattice is expected to immitate the presence of 
impurities in real physical systems. We usually appeal to the Harris crite- 
rion^ to answer the question of how such randomisation affects the critical 
exponents.^ If a > in the pure system, quenched disorder is deemed 
relevant and the critical exponents change as the disorder is increased. On 
the other hand, if a < in the pure model, then disorder of this type does 
not change the critical behaviour and the exponents are unaltered. In the 
case of a = 0, which as we have seen describes the Ising model in two 
dimensions, the Harris criterion does not provide a clear answer. 

For this reason, the bond- and site-diluted Ising models in two dimen- 
sions have been controversial over the years. The notion that some of the 
leading critical exponents change as the lattice structure is randomised, but 
where combinations which appear in terms of the correlation length such 
as j3/v and 7/1/ are unchanged, became known as the weak universality 
hypothesis. 

Gradually this gave way to the strong universality hypothesis which is 
now mostly believed to hold, although agreement is not universal.^ This 
predicts that the two-dimensional diluted models have the same leading 
critical exponents as in the pure case, but that there are multiplicative 
logarithmic correction s] 47 ! 83 ! 

a = 0, $ = ~^, 7=g, * = 0, v=\, 77 = 0. (1.193) 

With q = 0,^ these correction exponents obey the scaling relations for log- 
arithmic corrections (|1.39|) - (|1.42j) . The remaining scaling relations (|1.43jl - 
(TP61) predict 

A = - — , £=--, a c =-l, J> c = 0. (1.194) 
16' 2 ' y ' 

Moreover Eq. (|1.141[) means that there is a double logarithm in the specific 
heaflSMD 

Coo(t,0) ~ In I lnt|. (1.195) 
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The power of the scaling relations for logarithmic corrections is well illus- 
trated in this model as Eq. (|1.39p connects the hitherto most elusive and 
controversial quantity d directly to other exponents, which are more clearly 
established. For a review of the two-dimensional disordered Ising model, 
see Ref . [53] . 

1.6.9. Ashkin-Teller Model in 2 Dimensions 

The new relation (|1.127[) also holds in the /V-colour Ashkin-Teller model 
which consists of N coupled Ising models, the leading exponents are the 
same as Eq. (|1.192|) and^MI 

N - n-l . 7(N-1) * . N-l . 

r> P = oT' 7 = 77T7 ST, = 0, 1/ = — -, 77 = 0. 



N — 2 8(n-2)' ' 4(iV-2)' N - 2 

(1.196) 

If <7 = 0, these values also support the scaling relations (|1.39[) - (|1.42[) . 
Eqs. (|1.43|) - (|1.46|) further lead to the predictions 

i 15 ?i — 1 n — 1 re — 1 

A=- T r, e= s, «c = -2 -, ^ = 0. 1.197 

8 n — 2 n — 2 n — 2 

1.6.10. Spin Models on Networks 

Here we consider another set of phase transitions which exhibit multiplica- 
tive logarithmic corrections to scaling, namely spin systems on scale free 
networks. There are both academic and practical motivations for studying 
critical phenomena on complex networks.^ Phenomena such as opinion 
formation in social networks^! are expected to be modelled by such sys- 
tems, but one is also interested in realistic physics on complex geometries, 
such as for integrated nanoparticle systems.^ Some complex networks are 
characterised by so-called scale-free behavior. The degree of a node in the 
network is the number of links emanating from it and is denoted by k. 
With P(k) symbolising the degree probability distribution (the likelihood 
that an arbitrary chosen node has a certain degree value), this is power-law 
in scale-free systems, 

P(k) - k~ x . (1.198) 

For the critical phenomena previously described, dimensionality and 
length scales play crucial roles. Near the critical points itself, order- 
parameter fluctuations tend to be strongly correlated, the correlation length 
diverges and the pair-correlation function changes from an exponential to 
a power-law. The FSS hypothesis (|1 . 1 14|) depends upon the ratio of length 
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scales and the dimensionality enters the scaling relations through hyper- 
scaling (fTTT5|) and (fL39]) . 

However this notion of metrics is lost when we move from a lattice 
substrate to a complex network. For Euclidean lattices the coordination 
number is dimension dependent (it is 2d for a c?-dimcnsional hypcrcube). 
For networks the coordination numbers become the degrees associated with 
nodes and A plays the role of d. The node degree distribution function ex- 
ponent A in Eq. (|1.198p controls the non- homogeneity manifest in a network 
due to its internal structure. This is a principal difference between the 
origin of logarithmic corrections on regular lattices and on networks. It 
turns out that if A exceeds a critical value A c = 5 the phase transition 
has the usual mean-field critical exponents. As A decreases the node-degree 
distribution becomes increasingly fat-tailed and the relative amount of high- 
degree nodes (so-called hubs) increases. This leads to non-trivial critical 
behavior. As a result, systems with intermediate degree distributions where 
A s < A < A c have critical exponents which are generally A-dependent. Here 
A s = 3. Decreasing A still further to A < A s , the system becomes ordered at 
any finite temperature. Here only an infinite temperature field is capable 
of destroying the order. 

Just as in the regular-lattice case, where logarithmic corrections arise 
at d = d c , so too in the network case can they emerge at a marginal value 
A = A c for a number of classical spin models. In Ref. [29] . the field and 
temperature dependencies of critical thermodynamic functions were inves- 
tigated for a system with two coupled order parameters on a scale-free 
network. Models of this type are frequently used to describe systems with 
two different types of ordering. Physical examples of such systems include 
ferromagnetic and antiferromagnetic, ferroelectric and ferromagnetic, struc- 
tural and magnetic ordering. A sociophysics application^ may be opinion 
formation where there is a coupling between the preferences for a candidate 
and a party in an election, for example. 

A comprehensive description the system with two coupled scalar order 
parameters on a scale-free network was given in Ref. [29], where Eq. (|1.44j> 
was derived. The leading exponents (there are no exponents associated 
with correlation functions or the correlation length) when A = 5 are^l 



a = 0, (3 = 1/2, 7 = 1, 6 = 3, e=- 



a c = 
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In addition, the logarithmic-correction exponents are^l 




-, d c = -1, A 



1 



2' 

(1.200) 



This completes our list of systems exhibiting multiplicative logarith- 
mic corrections to scaling. The list is not exhaustive and the reader is 
invited to test the validity of the logarithmic scaling relations in other 
models which exhibit these phenomena. For example, three-dimensional 
anisotropic dipolar ferromagnets, in which spatial and spin degrees of free- 
dom are linked, are experimentally accessible syetems which exhibit log- 
arithmic corrections^ Systems with tricriticality also involve important 
logarithmic factors in three dimensions 

The four-dimensional diluted Ising model offers an example of a system 
with logarithmic corrections to leading scaling behaviour which is not power 
law. The XY model in two dimensions is an analagous example for a system 
with an infinite-order phase transition. 

1.7. Conclusions 

Over the past few years, a set of relations which link the exponents of log- 
arithmic corrections to scaling at higher-order phase transitions has been 
developed. These are the logarithmic analogues of the famous scaling re- 
lations between the leading critical exponents which were developed in the 
1960's and which are now pivotal in modern statistical mechanics as well 
as in related areas such as lattice quantum field theory. In this review, the 
logarithmic scaling relations arc presented and tested in a number of mod- 
els. In cases where there are gaps in the literature these scaling relations 
allow us to make predictions. 

With hindsight it is perhaps surprising that these logarithmic scaling 
relations have not been developed earlier. As mentioned in the text, while 
the relations between static exponents may be derived from a suitably mod- 
ified Widom hypothesis, this is not the case for the remaining exponents. 
One case which may have hindered earlier attempts to develop logarithmic 
scaling relations is that of the Ising model in two dimensions. There, only 
the specific-heat dependency on the reduced temperature has a logarithm. 
The fact that there appears, at first sight, to be no other logarithms to 
relate the specific-heat logarithm to, may be a reason why scaling relations 
for logarithmic corrections have not been developed earlier. 

Another stymying factor for an earlier development of the logarithmic 
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scaling relations is the rather enigmatic, new, logarithmic-correction critical 
exponent q which characterises the FSS of the correlation length. It was 
previously thoughtPSl that the finite-size correlation length could not exceed 
the actual length of the system. However, Eq. (|1.39j> shows not just that 
q may exceed zero, but that it is universal. This is clear because the 
exponents a and i>, with which q is related, are universal. This means that 
(a) sensible definitions of the finite-size correlation length and (b) sensible 
boundary conditions must respect the universal value of q. Any definition 
of the correlation length or any implementation of boundary conditions not 
respecting the universality of q can only be, in some sense, artificial. 

At the upper critical dimension, where q is non-trivial, the leading crit- 
ical exponents take on their mean-field values. Therefore they cannot be 
used to distinguish the universality class (e.g., as we have seen, the leading 
critical exponents for the 0(n) model in four dimensions are all the same). 
Instead, the exponents of the multiplicative logarithmic corrections may be 
used, so that these have a similar status there to the leading exponents for 
d < d c . 

The logarithmic-correction exponents given and derived above are listed 
in Tables 11.11 and 11.21 In the Table 11.21 bold-face symbols indicate the 
current gaps in the literature, which are filled by the theory presented 
herein and which now require independent verification. 
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Appendix A. Homogeneous Functions 

A function of a single variable f(x) is said to be homogeneous if multiplica- 
tive rescaling of x by an amount A results in a multiplicative rescaling of 



March 4, 2013 13:42 



World Scientific Review Volume - 9in x 6in 



RKchapter 



42 R. Kenna 



Tabic 1.1. Leading critical exponents for the models discussed herein. 
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f(x) by a factor <?(A), where g is a function of A only: 

f(Xx)=g(X)f(x). (A.l) 

Thus power-laws are homogeneous functions, but logarithmic functions are 
not. Indeed, if (|A.1[) holds, it turns out that both f{x) and g(x) have to be 



Table 1.2. Logarithmic-correction critical exponents for a variety of models. Bold face indicates gaps in the literature filled by the theory 
presented herein - i.e., bold- face values are predictions coming from the scaling relations for logarithmic corrections and remain to be verified 
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power laws. To see this, consider a second rescaling - this time by a factor 
\i. Rescaling first by A and then by /i leads to 

f(fiXx) = g(ti)f(\x) = g(fi)g(X)f(x), (A.2) 

while rescaling by the combined factor (jlX yields 

f(fiXx) = g(jiX)f(x). (A.3) 

Equating the right-hand sides of these two equations gives 

g(fiX) = g(tx)g(X). (A.4) 

Differentiating this with respect to fi gives Xg'(/j,X) = g'(fi)g(X), which, 
after chosing fi = 1 gives 

Xg'(X)=pg(X), (A.5) 

in which p = g'(l). Solving for g then gives g(X) oc X p , so that g(X) is 
power -law. Inserting this into Eq. (|A.lj) and chosing A = x -1 then gives 
/(l) = x~Pf(x), or 

f{x) ex x p , (A.6) 

so that f(x) is also power-law. 

A function of two variables f(x,y) is called homogeneous if 

f(Xx,Xy)=g(X)f(x,y). (A.7) 

Following similar considerations one can easily show that g(X) is again 
power-law. A more generalised homogeneous function obeys 

f(X r x,X s y) = XPf(x,y). (A.8) 

Renaming X p — > X one arrives at the following property for a generalised 
homogeneous function: 

f(X a x,X b y) = Xf(x,y). (A.9) 
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